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SUMMARY 


In this reporting period, further experimentations were made to 
improve the design and fabrication techniques of the integrated sole. 

The sole design is shown to be related to the foot position 
requirements and the actual shape of the foot including presence of 
neurotropic ulcers or other infections. Factors for consideration were, 
heel pitch, balance line, and rigidity conditions of the foot. Machining 
considerations were also part of the design problem. Among these 
considerations, widths of each contour, tool motion, tool feed rate, depths 
of cut, and slopes of cut at the boundary were the key elements. 

The essential fabrication techniques evolved around the idea of 
machining a mold then, using quick-firm latex material, casting the sole 
through the mold. Two main mold materials were experimented with: plaster 

and wood. Plaster was very easy to machine and shape but could barely 
support the pressure in the hydraulic press required by the casting 
process. Wood was found to be quite effective in terms of relative cost, 
strength, and surface smoothness except for the problem of cutting against 
the fibers which could generate ragged surfaces. 

This report also discusses our programming efforts to convert our 
original dBase programs into C programs so that they could be executed on 
the SUN Computer at North Carolina State University. 
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2 - DESIGN CONSIDERATIONS 


There are two major categories of design considerations: 

(1) Orthopedic /Therapeutic 

(2) Machining 

The Orthopedic/Therapeutic requirements arise from physiological 
disorders of the patient and dictate the types of sole needed to alleviate 
or prevent further damage to his foot. The Machining requirements are 
strictly technological requirements to achieve the design objectives of the 
sole. 


2.1 Orthopedic/Therapeutic Requirements 

People who may benefit from custom-molded shoes can be a healthy 
person with hard- to- fit feet or, as usually is the case, a person with 
severe physiological disorders of lower limbs. In the majority of cases, 
the custom shoe candidate is at least partially disabled. He or she may 
have been born with a birth defect such as cerebral palsy, or club feet. 
He or she may have had diseases such as rheumatoid arthritis, diabetes, or 
polio. Muscular dystrophy, multiple sclerosis, or charcot disease may have 
been his or her problem. Aging also can contribute to foot imbalance, 
bunions, or hammertoes. Regardless of etiology, this person cannot 
tolerate commercial footwear, and even extra-depth shoes with special 
inserts may only offer minimum reduction in discomfort and pain. It is 
this particular segment of population for whom custom molded shoes are a 
requirement that is addressed in this NASA supported project. 
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the length of the tibia but, when conditions of foot inflexibility and 
severe proration or supiration exist, it should be selected appropriately 
and used to orient the last before one can proceed with the design of the 
sole. 

Foot rigidity : Generally speaking if the foot is flexible then a 
certain amount of cast inversion or eversion is tolerated. This means that 
the angle of inversion or eversion could be specified by the user. If the 
foot is inflexible or deformation conditions such as cavus foot, 
equino-varus , or considerable tibia varum exist then the balance line is 
drawn while the foot is in its natural orientation. 

The discussion of heel pitch, balance line, and foot rigidity given 
above has been communicated to the research group at North Carolina State 
University for incorporation into their LASTMOD* design software. See 
Appendix 1. 

2.2 Machining Requirements 

The basic concept of making the sole for the shoe is to generate a 
mold cavity with its bottom surface being an exact reproduction of the 
plantar surface of the foot in an upside-down position. This concept is 
illustrated in figure 1. 

Plantar surface in upside down position 


a) machining of mold b) casting of sole using liquid 

latex 

Figure 1 - Concept of mold making 
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Custom molded shoes require a last which remains inside the shoes from 
the beginning to the end of the shoe-making process. The last is the most 
important element of the entire process and must be precisely designed and 
fabricated. The responsibility for last design and last machining is with 
the research group at North Carolina State University, while the design and 
manufacture of insoles are the responsibility of the University of Missouri- 
Columbia. 

As far as the design of the soles are concerned, the foot must first 
be scanned in its most correct position. The correct position is typically 
a "neutral" position in which there is minimum stress to the foot. 
Opinions vary between proponents of no-load position and proponents of 
weight-bearing position but the general practice is the seated position for 
which the foot is at rest. The knees should be flexed at 90° unless there 
is limited flexion in the knees. The position of the tibia is also an 
important consideration, as well as the orientation of the ankle joints 
because of their influence on the heel pitch. Reference [1] provides ample 
discussion on the factors just mentioned. 

In general, there are three basic design criteria for the soles: heel 

pitch, balance line, and rigidity condition of the foot. 

Heel pitch : It is specified according to aesthetic requirements (ie., 
male or female) and anatomical features such as leg-length discrepancy. 
There exist general heel pitch values (1/2" to 1 3/4"), and they should be 
used as much as possible. 

Balance line : The balance line can be drawn on the last to indicate 

the amount of correction desired in the shoe. It is usually vertical along 
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Before a milling machine is used to mill the mold cavity, a number of 
design problems must first be resolved. They are listed as follows: 
undercut 
high edges 
cutter interference 

2*2.1 Undercut 

This problem is illustrated in figure 2. To generate the bottom 
surface for the mold cavity, a tool must be programmed to follow each 
contour of the foot. However, if we allow the tool to go below the maximum 
width (shown as dash line in figure 2) then some material may be removed 
from the foot body, thus distorting its final shape. 



I 


to leg 


Figure 2 - Foot contour and problem of undercut 
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The solution, therefore, is to program the tool motion so that it will 
never go below the points of maximum width as shown in figure 2. The 
implication of this strategy is a cut surface through the foot which is 
neither planar nor constant width as shown in figure 3. The fact that this 
surface is everywhere as wide as the foot can be means that the foot will 
not feel squeezed inside the shoe and, therefore, the fit of the shoe would 
be enhanced. 



line of maximum width 


Figure 3 - Line of maximum width 


2.2.2 High edges 

When the points of maximum width happens to be well above the bottom 
surface of the foot, high edges may result as shown in figure 4. 
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High edges present two problems. First, they require very thin 
cutters which may break easily during the machining phase. Second, the 
casting latex may not be able to reach the bottom of these edges, thus 
creating a ragged rather than smooth edge. 

These problems have been avoided in our fabrication experiments by 
arbitrarily limiting the slope to no more than 30 degrees. 

2.2.3 Cutter Interference 

Cutter interference is defined as a problem of inadvertently removing 
some portions of the part surface as due to the size of the cutter. This 
problem is illustrated in figure 5. 



Figure 5 - Cutter Interference 

Two types of cutter interference have been identified: 
circumferential and longitudinal. 
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Circumferential Interference 


This problem occurs when the tool is programmed to move along each 

contour of the shoe last as shown in figure 5. 

The approach was to calculate the amount of tool lift necessary to 

avoid the interference, very much along the line used by Saunders and 

Vickers [2] except for the fact that the tool is always vertical while, in 

their case, the tool is radially oriented. The tool lift was calculated as 

follows. As shown in figure 6, let's consider a point of interest, say t^, 

and its immediate neighbors t. +1> and t^ respectively before and behind 

it. The tool is assumed to move from t^_^ to to t^ + ^. Interference is 

detected if the outline of the tool centered at t, intersects with segment 

t.t..,. The maximum vertical difference between segment t.t. and the 
l l+l 1 1 “*' 1 

tool contour (a circle) is taken as the necessary tool uplift required to 
avoid circumferential interference. 



Figure 6 - Tool lift to avoid circumferential interference 
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Longitudinal Interference 


This problem occurs when the tool is programmed to move along the 
length of the object. In our case the tool spends most of its time cutting 
across the object (ie., circumferential cuts), and very little time along 
the length of the object except for advancing the tool forward for the next 
cut. Because this motion corresponds with the conventional feed rate, the 
problem of longitudinal interference reduces to a problem of surface 
roughness left between two consecutive tool passes. Because of the nature 
of the cut, longitudinal interference can be ignored without much effect on 
the surface to be generated. 

3 - MACHINING WORK 


To implement the design objectives discussed in section 2, a series of 
computer programs written mostly in dBase language were developed as shown 
in figure 7. figures 8 to 11 show respectively the 4-axis CNC milling 
machine used to machine the mold cavities, a mold in plaster, a mold in 


wood, and samples of casted insoles. 



Select data 
points below 
max width 


Increase 
resolution by 
in t er p o la t ion 


Continued on next page 
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DFLATED1.PRG 



insert another contour 
DINTER2.PRG between every two consecutive 
' [ contours 


Reduce the high edges 


DS0LET02 •DBFi 


DINTER3.PRG 


insert four more contours 
between 10.00 and 10.25 


DS0LET03 .DBF] 


DREB01.PRG 


Form the zig-zag cutting 
sequence 


DS0LET04 .DBF) 


[DS0LET06.DBF 




DINFERI.PRG 




Tool interference 


DS0LET07 .DBF 




DNCT07 .PRG 




Convert DS0LET07.DBF 
to NC Program 


DNCT07.PRG 


DI0S1.PRG 


Add heel height to sole 
surface . 


DS0LET05.DBF 



Data File 



Program 


Figure 7 - dBase programs to generate machine codes 
for machining mold cavity 










Figure 8 4 -Axis CNC Milling Machine Used in 

ORIGINAL PAGE IS the Production of the Mold Cavity. 

OF POOR QUALITY 



Figure 9 An example of Plaster Mold 


ORIGINAL page* 

BLACK AND WHITE PHOTOGRAPH 
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Figure 10 An Example of Wood Mold 

ORIGINAL PAGE IS 
OF POOR QUALITY 



ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


Figure 11 Sample Insoles 
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4 - SOFTWARE CONVERSION 


Since the final product of this project is a custom footwear design 
system based on the SUN computer, all software developed by the various 
parties must be converted to C programs. 

Figure 12 shows a sequence of three C programs designed to bring the 
work done by the group at the University of Missouri-Columbia into the 
LASTMOD' program at North Carolina State University. 



Machining Codes for 
milling machine 



Data File 


Program 


Figure 12 - Sequence of C programs to make the current work 
compatible with the SUN computer 
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DETACH. C 


This program searches for the maximum width of each foot contour, 
discards all data points above this width, reduces the high edges to an 
arbitrary 20-degree slope, and rearranges the data points so that a 
continuous zig-zag tool motion can be obtained. 

INSOLE. C 

This program finds the highest points at the heel and at the ball of 
the foot, adds user-specified heights to these points, for example 1 3/4 

inch to heel and 1/2 inch to ball, then determines the thicknesses of the 
sole at all data points. 

NCCODE.C 

This program generates the machining codes for the milling machine. 
For this, it must add G, M, S, and F codes to each motion of the tool. 


The above programs have been installed in the LASTMOD software at the 
last group meeting at North Carolina State University on May 23, 1990. 
Their listings are shown in appendix 3. 
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APPENDIX 1 


UNIVERSITY OF MISSOURI-COLUMBIA 


Department of Industrial Engineering 

113 Electrical Engineering Building 
Columbia, Missouri 65211 
Telephone (314) 862-2691 


MEMORANDUM 


TO 

FROM 

SUBJECT 

DATE 


Dr. D. McAllister fj»/7 
Han P . Bao \jff // 

Sole Design : features to be included in LASTMOD 
June 19 , 1989 


(cc : Drs Farmer, Luo, Rasdorf, Sanii ) 


This memorandum contains three parts : 

1- fundamentals of sole design 

2- features to be included in LASTMOD 

3- types of output files to be given to UMC 

Part I outlines the essential specifications of last position 
before the design of the matching sole can be pursued. 

Part II indicates the desirable features ( ie. new features ) 
to be included in LASTMOD as part of the sole design process. 

Part III specifies the types and formats of the output files 
that your research group should provide us. 


PART I 

For the sole to be designed properly, we require two pieces of 
information : 

1- heel pitch 

2- balance line 


Figure 1 depicts these two features 



fit civ 



Line. 


an equal opportunity institution 


Pt A : lowest pt of 
heel 

Pt B : lowest pt of 
ball of foot 



4- „Once the heel pitch is selected then redisplay the 

foot with the selected heel pitch for visual check. ^Ef 
the view is satisfactory, then 1 done 1 option is - 
selected, and the panel in step 2 is shown again. 


5- If 1 balance line 1 option in step 2 is selected then 
display the following screen : 



yes 

no 

flexibility 


□ 

inversion angle 

(default is 


zero) 

eversion angle 

Tdefault is 


zero) 

done j 




When 1 done 1 is picked in step 5, the screen in step 
2 should be shown again. 


In summary, the four basic pieces of information required for the sole 
design process are : 

1- heel pitch 

2- inversion angle 

3- eversion angle 

4- foot flexibility ( yes or no ) 

They will be used to determine a reference plane below which the foot 
would be shown and used in the derivation of the mold geometry for 
casting the sole. 



The heel pitch is specified according to aesthetic requirement ( male 
or female fashion ) and anatomical features such as leg-length dis- 
crepancy ( LLD ). There exist general heel pitch values, and they 
should be indicated on the computer screen for selection. 

The balance line is a line drawn on the back of the last for purpose 
of balancing the cast after the addition of the heel. This balance 
line is usually vertical along the length of the tibias but, when 
conditions of foot inflexibility and severe pronation or supination 
exist, it should be selected appropriately and used to orient the last 
before one can proceed with the design of the sole. 

In general, if the foot is flexible then a certain amount of cast 
inversion or eversion is tolerated. This means that the angle of 
inversion or eversion must be specified interactively by ther user. If 
the foot is inflexible or deformation conditions such as cavus foot, 
equino-varus , or considerable tibia varum exist then the balance line 
is drawn while the foot is in its natural orientation. 


PART II 


The additional features for LASTMOD will be : 

1- in MAIN menu, add option 1 sole design * 

2- Design a new panel labelled 'sole design cycle^ to 

look like this. 



note : this panel should appear after the 'sole design' 
option in the main menu is selected. 

3- If ' heel pitch section 1 in step 2 is selected then display 
medial view of foot and sub-menu as shown below. 




lowest pt of 
heel 

lowest pt of 
ball of foot 



flYehs 


heel 

pitch 

1/2 

- 3/4 

L 

1 1/2 
2 

select 

done 


Once the heel pitch is selected then redisplay the 
foot with the selected heel pitch for visual check, if 
the view is satisfactory, then 'done ' option is '■ 
selected, and the panel in step 2 is shown again. 

If ' balance line * option in step 2 is selected then 
display the following screen : 



yes 

no 

flexibility 


□ 

inversion angle 


(default is 

eversion angle 


zero) 
Tdefault is 

done 


zero) 


When 1 done ' is picked in step 5, the screen in step 
2 should be shown again. 


In summary, the four basic pieces of information required for the sole 
design process are : . 

1- heel pitch 

2- inversion angle 

3- eversion angle 

4- foot flexibility ( yes or no ) 

They will be used to determine a reference plane below which the foot 
would be shown and used in the derivation of the mold geometry for 
casting the sole. 




PART III 


Two types of files are required : 

a- parameter file to contain last ID number, heel pitch, 
inversion angle, eversion angle, and foot flexibility. 

b- Points coordinate file to contain the geometry of the 
lower portion of the foot ( see examples below ) 


Example 1 : normal foot, inversion angle = eversion angle = JD 



Example 2 : pronated foot, flexible, eversion angle=5 deg 
inversion = 0 



Examp le 3 : supinated foot, flexible, inversion angle = 5 
deg., eversion angle = 0 deg. 



inversion angle 


******************* 
close databases 
select 2 
use side 
select 1 
use foot2 
rx=x 


PROGRAM NAME : Sl.PRG 




» 


maax=x 

maay=y 

maaz=z 

miix=x 

miiy=y 

miiz=z 

skip 

do while .t. 
select 1 
if .not. (eof () ) 
if xorx+0.5 
if y>maay 
maay=y 
maax=x 
maaz=z 
skip 
loop 
else 

if ycmiiy 
miix=x 
miiy=y 
miiz=z 
skip 
loop 
else 
skip 
loop 
endif 
end if 
else 

select 2 
append blank 

replace lax with maax,lay with maay,laz with maaz 
replace smx with miix,smy with miiy,smz with miiz 
select 1 
rx=x 
maax=x 
maay=y 
maaz=z 
miix=x 
miiy=y 
miiz=z 
loop 
endif 
else 

select 2 
append blank 

replace lax with maax,lay with maay,laz with maaz 
replace smx with miix,smy with miiy,smz with miiz 
exit 
endif 
enddo 


******************* PROGRAM NAME : SOLE2.PRG ************************* 
close databases 
select 1 

use foot2 * 

rx=x 

select 2 

use side 

select 3 

use sole 

f lagl=0 

f lag2=0 

f lag3=0 

do while .t. 

select 1 

if x<>rx+0 . 5 . and. ( .not.eof () ) 
if y=side->lay 
f lagl=l 
endif 

if y=side->smy 
f lag2=l 
flag3=l 
endif 
do case 

case f lagl=0 . and . f lag2=0 . and . f lag3=0 
select 1 
skip 
loop 

case flagl=l.and. flag2=0.and. flag3=0 
select 3 
append blank 

replace sox with foot2->x,soy with foot2->y,soz with foot2-> 

select 1 

skip 

loop 

case flagl=l.and. flag2=l.and. flag3=l 
select 3 
append blank 

replace sox with foot2->x,soy with foot2->y,soz with foot2-> 

flag3=0 

select 1 

skip 

loop 

case f lagl=l . and . f lag2=l . and . f lag3=0 
select 1 
skip 
loop 


endcase 

else 

if . not . (eof ( ) ) 
flagl=0 
flag2=0 
rx=x 

select 2 
skip 
loop 
else 
exit 


endif 
endif . 
enddo 


******************* PROGRAM NAME : DINTER1 . PRG ************************* 

* PURPOSE: TO INSERT FOUR MORE SLICES IN BETWEEN EVERY TWO CURRENT * 

* SLICES BY INTERPOLATION ( 0.5->0.1 INCH ) * 

* PROGRAMMER: PAUL N. LIN DATE: DEC. 6 1989 * 

*****************************************************i******************** 

close databases 
select 1 
use dsole2 
select 2 
use dsole 


rx=sox 

inx=sox 


nol=recno() 
do while .t. 

count next 25 for sox=rx to ccl 
goto nol 

count next 50 for sox=rx+0.5 to cc2 
do dsubintl with nol,ccl,cc2 
if soxoiO.O 
rx=rx+0 . 5 
loop 
else 

do while .t. 

if .not.eof() 
select 1 
append blank 

replace sox with dsole— >sox, soy with dsole— >soy, soz with ; 
dsole— >soz 

select 2 

skip 

loop 

enddo 

exit 

endif 

enddo 

return 


******************* PROGRAM NAME : DSUBINT1. PRG ************************* 

* PURPOSE: THE SUBROUTINE OF DINTER1 . PRG TO ADD SLICES IN BETWEEN * 

* CURRENT SLICES ( 0.5 -> 0.1 INCH ) * 

* PROGRAMMER: PAUL N. LIN DATE: ‘DEC. 6, 19 89 * 

*******************************************************************^^^^^^^^ 
parameter 001,001,002 

nno=nol 

head=max ( ccl , cc2 ) 
tail=min (ccl , cc2 ) 
cut=mod (head , tail ) 
if cut=0 
pa=0 
else 

if head=ccl 
pa=l 
else 

if head=cc2 
pa=2 
endif 
endif 
endif 
do case 

case pa=0 
num=0 

do while num<5 
go nol 
ppl=ccl 

irx=sox+ ( 0 . l*num) 
do while ppl>0 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+(ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl-l 
select 2 
skip 
loop 
enddo 
num=num+l 
loop 
enddo 

nol=recno ( ) 
return 
case pa=l 
go nol 

odd=mod(ccl,2) 
do case 

case odd=0 

midd=(ccl/2) 


rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsole->sox, soy with dsole->soy , soz with 
dsole->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
PPl-1 

irx=sox+ (0 . l*num) 
do while ppl<=ccl 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
else 

noo=recno ( ) 

oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz ) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 

ppl=ppl+l 

select 2 
skip 
loop 
endif 
enddo 
num=mm+l 


loop 

enddo 

nol=recno ( ) 
return 
case odd=l 

midd= (int (ccl/2) ) 
rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsole->sox, soy with dsole->soy , soz with 
dsole->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
'ppl=l 

irx=sox+ ( 0 . l*num) 
do while ppl<=ccl 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
else 

noo=recno ( ) 

oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 



select 2 
skip 
loop 
endif 
enddo 
num=num+l 
loop 
enddo 

nol=recno ( ) 
return 

endcase 
case pa=2 
go nol 

odd=mod ( ccl , 2 ) 
do case 

case odd=0 

midd=(ccl/2) 

rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsole->sox, soy with dsole->soy, soz with 
dsole->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
ppl=l 

irx=sox+ ( 0 . l*num) 
do while ppl<=cc2 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly~(dy-oy)/5 
lz= (dz-oz ) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
else 

if ppl=midd+l 
select 2 
skip -1 
endif 

noo=recno ( ) 


oy=soy 
oz=soz 
sk=noo+cc2 
go sk 
dy=soy 
dz=soz ! 
ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
endif 
enddo 
num=num+l 
loop 
enddo 

nol=recno() 
return 
case odd=l 

midd= ( int ( ccl/2 ) ) +1 
rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsole->sox, soy with dsole->soy , soz with 
dsole->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
ppl=l 

irx=sox+ ( 0 . l*num) 
do while ppl<=cc2 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz ) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 



ppl=ppl+l 
select 2 
skip 

loop * 

else 

if ppl=midd+l 
select 2 
skip -1 
endif 

noo=recno ( ) 

oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly *num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
endif 
enddo 
num=num+l 
loop 
enddo 

nol=recno ( ) 
return 

endcase 

endcase 



******************* PROGRAM NAME : DINTER2 . PRG ************************ 

* THIS PROGRAM IS USED TO INSERT ONE MORE LINE IN BETWEEN EVERY TWO * 

* SLICES AND MAKE THE DISTANCE BETWEEN EVERY TWO CONTURE TO BE 0.05 * 

* INCH. ( DSOLE . DBF -> DS0LET02 . DBF ) ‘ * 

* PROGRAMMER : PAUL N. LIN DATE : DEC. 6, 1989 * 

***********************************, 1 .^^^^^^.^**^^^.^**^^^^^^^^^^.^^^^^^^^^^^^ 

close databases 
select 1 
use dsoleto2 
select 2 
use dsole2 
rx=sox 
inx=sox 
nol=recno() 
do while .t. 

count next 25 for sox=rx to ccl 
goto nol 

count next 50 for sox=rx+0.1 to cc2 
do dsubint2 with nol ,ccl,cc2 
if soxoiO.O 
rx=rx+0 . 1 
loop 
else 

do while .t. 

if .not.eof() 
select 1 
append blank 

replace sox with dsole2— >sox, soy with dsole2->soy , soz with ; 
dsole2->soz 

select 2 

skip 

loop 

enddo 

exit 

endif 

enddo 

return 



******************* PROGRAM NAME : DSUBINT2 . PRG ************************ 

* USAGE : A SUBROUTINE OF DINTER2 . PRG TO MAKE THE DISTANCE BETWEEN EVERY * 

* TWO SLICES TO BE 0.05 INCH. * 

* PROGRAMMER : PAUL N. LIN DATE : DEC.* 6, 19 8 9 * 

A*************************************************^^^^^^^^^*^*^ 

parameter nol,ccl,cc2 
nno=nol 


head=max(ccl, cc 2 ) 
tail=min (ccl , cc2 ) 
cut=mod (head, tail) 
if cut=0 


pa=0 

else 

if head=ccl 
pa=l 
else 

if head=cc2 
pa=2 
endif 
endif ' 
endif 
do case 

case pa=0 
num=0 

do while num<2 
go nol 
ppl=ccl 

irx=sox+ (0. 05*num) 
do while ppl>0 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/2 
lz=(dz-oz )/2 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 

ppl=ppl-l 
select 2 
skip 
loop 
enddo 
num=num+l 
loop 
enddo 

nol=recno ( ) 
return 
case pa=l 
go nol 

odd=mod (ccl , 2 ) 
do case 

case odd=0 

midd=(ccl/ 2 ) 



rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsole->sox, soy with dsole->soy , soz with 
dsole->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
ppl=l 

irx=sox+ (0 . l*num) 
do while ppl<=ccl 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy - 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
else 

noo=recno ( ) 

oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly-(dy-oy) /5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*nuin) 
select 1 
append blank 

replace sox with irx ; soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
endif 
enddo 
num=num+ 1 



loop 

enddo 

nol=recno () 
return 
case odd=l 

midd=(int (ccl/2) ) 
rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsole->sox, soy with dsole->soy , soz with 
dsole->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
ppl=l 

irx=sox+ ( 0 . l*num) 
do while ppl<=ccl 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz ) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
else 

noo=recno ( ) 

oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz ) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+I 



I 


I 

I 

I 

I 

I 

a 
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select 2 
skip 
loop 
endif 
enddo 
num=num+ 1 
loop 
enddo 

nol=recno ( ) 
return 

endcase 
case pa=2 
go nol 

odd=mod ( ccl , 2 ) 
do case 

case odd=0 

midd= (ccl/2 ) 
rr=ccl 

do while rr>0 
select 1 
. ' append blank 

replace sox with dsole->sox, soy with dsole->soy , soz with 
dsole->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
ppl=l 

irx=sox+ ( 0 . l*nuxn) 
do while ppl<=cc2 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz ) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
else 

if ppl=midd+l 
select 2 
skip -1 
endif 

noo=recno ( ) 



oy=soy 
oz=soz ' 
sk=noo+cc2 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
endif 
enddo 
num=num+l 
loop 
enddo 

nol=recno ( ) 
return 
case odd=l 

midd=(int(ccl/2) )+l 
rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsole— >sox, soy with dsole— >soy , soz with ; 
dsole->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
ppl=l 

irx=sox+ ( o . l*num) 
do while ppl<=cc2 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz ) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 



ppl=ppl+l 
select 2 
skip 

loop * 

else 

if ppl=midd+l 
select 2 
skip -1 
endif 

noo=recno ( ) 

oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly=(dy-oy)/2 

lz= (dz-oz) /2 

go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
endif 
enddo 
mxm=num+l 
loop 
enddo 

nol=recno () 
return 

endcase 

endcase 


******************* PROGRAM NAME : DINTER3 . PRG ************************* 

* PURPOSE: TO INSERT FOUR MORE SLICES IN BETWEEN 10.00 and 10.25 * 

* ( 0 . 2 5->0 . 05 INCH ) * 

* PROGRAMMER: PAUL N. LIN DATE: DEC. 6, 1989 * 

******************************************************* a ***** a ** +a ^^^ + ^ a ^ + 
close databases 

select 1 
| use dsoleto3 
select 2 
use dsoleto2 
i do while .t. 
hx=sox 
if hxoiO.OO 
select 1 
append blank 

| replace sox with dsoleto2— >sox, soy with dsoleto2->soy , soz with ; 

dsoleto2->soz 

select 2 
skip 
loop 
end if 
exit 
enddo 
rx=sox 
inx=sox 
nol=recno ( ) 

count next 25 for sox=rx to ccl 
goto nol 

count rest for sox=rx+0.25 to cc2 
do dsubint3 with nol,ccl,cc2 
do while .t. 

if .not.eof() 
select 1 
append blank 

replace sox with dsoleto2— >sox, soy with dsoleto2— >soy , soz with ; 
dsoleto2->soz 

select 2 
skip 
loop 
end if 
exit 
enddo 

close databases 
return 


******************* PROGRAM NAME : DSUBINT3 . PRG **************** 

* PURPOSE: THE SUBROUTINE OF DINTER3 . PRG TO ADD FOUR MORE SLICES 

* IN BETWEEN 10.00 AND 10.25 ( 0.25 -> 0.0.05 INCH ) 

* PROGRAMMER: PAUL N. LIN DATE: DEC. 6, 1989 

a***********************************************^;^;^^^^^^ 

parameter nol,ccl,cc2 


********* 

* 

* 

* 

********* 


nno=nol 


head=max (ccl , cc2 ) 
tail=min ( ccl , cc2 ) 
cut=mod (head , ta il ) 
if cut=0 


pa=0 

else 

if head=ccl 
pa=l 
else 

if head=cc2 
pa=2 
end if 
endif - 
endif 
do case 

case pa=0 
num=0 

do while num<5 
go nol 
ppl=ccl 

irx=sox+ ( 0 . l*num) 
do while ppl>o 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz=(dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 

ppl=ppl-l 

select 2 
skip 
loop 
enddo 
num=num+l 
loop 
enddo 

nol=recno ( ) 
return 
case pa=l 
go nol 

odd=mod(ccl, 2) 
do case 

case odd=0 

midd=(ccl/2) 



rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsoleto2->sox, soy with dsoleto2->soy , 
soz with dsoleto2->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
ppl=l 

irx=sox+ ( 0 . l*num) 
do while ppl<=ccl 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly~ (dy-oy) /5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
else 

noo=recno ( ) 

oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly= (dy-oy) /5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly *num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
endif 
enddo 
num=nuirt+l 



case 


loop 

enddo 

nol=recno ( ) 

return 

odd=l 

xnidd= ( int (ccl/2 ) ) 
rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsoleto2->sox, soy with dsoleto2->soy , 
soz with dsoleto2->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
nuxn=l 

do while num<5 
go nol 
-ppl=l 

irx=sox+ (0 . l*num) 
do while ppl<=ccl 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy ; 
dz=soz • 
ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
else 

noo=recno ( ) 

oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 



select 2 
skip 
loop 
endif 
enddo 
num=num+l 
loop 
enddo 

nol=recno ( ) 
return 

endcase 
case pa=2 
go nol 

odd=mod ( ccl , 2 ) 
do case 

case odd=0 

midd=(ccl/2) 

rr=ccl 

do while rr>0 
select 1 
"append blank 

replace sox with dsoleto2->sox, soy with dsoleto2->soy , 
soz with dsoleto2->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
ppl=l 

irx=sox+ ( 0 . l*num) 
do while ppl<=cc2 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 

ppl=ppl+l 

select 2 
skip 
loop 
else 

if ppl=midd+l 
select 2 
skip -1 
endif 

noo=recno ( ) 



oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz) /5 
go noo 

iry=soy+(ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 

ppl=ppl+l 

select 2 
skip 
loop 
endif 
enddo 
'num=num+l 
loop 
enddo 

nol=recno ( ) 
return 
case odd=l 

midd=(int (ccl/ 2 ) )+l 
rr=ccl 

do while rr>0 
select 1 
append blank 

replace sox with dsoleto2->sox, soy with dsoleto2->soy, 
soz with dsoleto2->soz 

rr=rr-l 
select 2 
skip 
loop 
enddo 
num=l 

do while num<5 
go nol 
ppl=l 

irx=sox+ ( o . l*num) 
do while ppl<=cc2 
if ppl<=midd 
noo=recno ( ) 
oy=soy 
oz=soz 
sk=noo+ccl 
go sk 
dy=soy 
dz=soz 

ly=(dy-oy)/5 
lz= (dz-oz ) /5 
go noo 

iry=soy+ ( ly*num) 
irz=soz+ ( lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 



ppl=ppl+l 

select 2 
skip 

loop * 

else 

if ppl=midd+i 
select 2 
skip -1 
endif 

noo=recno ( ) 

oy=soy 

oz=soz 

sk=noo+cc2 

go sk 

dy=soy 

dz=soz 

ly=(dy-oy)/5 

lz= (dz-oz) /5 

go noo 

iry=soy+ ( ly*num) 
irz=soz+ (lz*num) 
select 1 
append blank 

replace sox with irx,soy with iry,soz with irz 
ppl=ppl+l 
select 2 
skip 
loop 
endif 
enddo 
num=num+l 
loop 
enddo 

nol=recno ( ) 
return 

endcase 

endcase 



;v*:;r**r******* program name : RE boi.pr G **************************** 

* : E?ER?^wf CO^oSRflsToS^N^f SHAPE “ DISTAHCE BETWEEN 

* PROGRAMMER : NIENSUNG LIN DATE - OCT *5 i q n q 

close*databases************************ ************************ ************ 

select 2 

use soleto3 

select l 

use soleto2 

rx=sox 

sh=0 

cc=0 

do while .t. 

if sox=rx . and . ( . not . eo f ( ) ) 
select 2 
append blank 

s!lectV° X With soleto2 - >s °x.soy with soleto2->soy , soz with soleto2-: 

skip 

loop 

else 

if sox=rx+0 . 05 . and . ( . not . eof f ) ) 
cc=cc+l 


skip 
loop 
else 
if 


sox=rx +0 
rx=sox 
bb=cc 
skip -l 
do while 
if 


1 . and . ( . not . eof ( ) ) 


.t. 

cc<>0 

1 select 2 

append blank 

repiace so X vith^sole^-^ox.soy with soleto2->soy,soz with soleto 2 ->sc 

cc=cc-l 
skip -l 
loop 
else 

select 1 
dd=recno ( ) +bb+l 
goto dd 
exit 
endif 
enddo 
loop 
else 

if eof() 
exit 
endif 
endif 
endif 
endif 
enddo 


************************** PROGRAM NAME : IOS.PRG ********************** 

* USAGE : TO ADD A PROPER THICNESS TO THE SOLE AND TURN IT OVER FOR THE * 

* NECESSITY OF MACHINING WHERE THE FOOT DATA IS OF 0.05 INCH * 

* CONTOUR DISTANCE. * * 

* PROGRAMMER : NIENSUNG LIN DATE : OCT. 5 1989 * 

************************************************************************** 
close databases 

select 1 
use soleto5 
select 2 
use soleto4 
rx=sox 

inc=(2-l)/200 

ba=-2 

do while ,t. 

if soxorx+0 . 05 . and. ( . not . eof ( ) ) 
iiz=ba-soz 
select 1 
append blank 

replace sox with soleto4->sox, soy with soleto4->soy , soz with iiz 
select 2 
skip 
loop 
else 

if (sox=rx+0 . 05) .and. (.not. eof () ) 
rx=sox 
ba=ba+inc 
loop 
else 

if eof ( ) 
exit 
end if 
endif 
end if 
enddo 

close databases 
return 



********************** PROGRAM NAME : FLATED1 . PRG ************************ 

* USAGE : TO FLAT THE EDGE OF THE SOLE AND MAKE IT SMOOTHER * 

* PROGRAMMER : NIENSUNG LIN DATE : OCT. 6 1989 * 

**************************************************************************** 
close databases 

select 2 
use soleto5 
ix=sox 
do while .t. 

if sox= ( ix+O . 05) . and. ( . not . eof ( ) ) 
skip -1 
iy=soy 
skip -1 
do while .t. 

if iy<0.0001 

if soy>(iy+0.5) 
fx=sox 
fy=soy 
fz=soz 
skip 

do while .t. 

if soxo (ix+0 . 05) 
sx=sox 
sy=soy 
sz=soz 

m=(sz-fz)/ (sy-fy) 
if m>l 

ssz=f z+ (sy-fy) 
replace soz with ssz 
endif 
fx=sox 
fy=soy 
f z=soz 
skip 
loop 
else 
exit 
endif 
enddo 
ix=sox 
exit 
else 

skip -1 
loop 
endif 
endif 

if iy>-0 .0001 

if soy<(iy-0.5) 
fx=sox 
fy=soy 
f z=soz 
skip 

do while .t. 

if soxo (ix+0. 05) 
sx=sox 
sy=soy 
sz=soz 

m=(sz-fz)/ (sy-fy) 
if m<(-l) 

ssz=fz- (sy-fy) 


replace soz with 
end if 
fx=sox 
fy=soy 
f z=soz 
skip 
loop 
else 
exit 
endif 
enddo 
ix=sox 
exit 
else 

skip -1 
loop 
endif 
endif 
enddo 
loop 
else 

if ( .not. eof ( ) ) 
skip 
loop 
else 
skip -1 
iy=soy 
do while .t. 
if iy<0.0001 

if soy> ( iy+O . 5) 
fx=sox 
fy=soy 
fz=soz 
skip 

do while .t. 

if ( .not. eof () ) 
sx=sox 
sy=soy 
sz=soz 

m=(sz-fz)/ (sy-fy) 
if m>l 

ssz=fz+ (sy-fy) 
replace soz with 
endif 
fx=sox 
fy=soy 
f z=soz 
skip 
loop 
else 

exit 

endif 

enddo 

exit 

else 

skip -1 
loop 
endif 
endif 

if iy>-0.0001 


ssz 


ssz 



if soy<(iy-0.5) 
fx=sox 
fy=soy 
fz=soz 
skip 

do while .t. 

if ( .not. eof ( ) ) 
sx=sox 
sy=soy 
sz=soz 

m= (sz-fz) / (sy-fy) 
if m<(-l) 

ssz=fz- (sy-fy) 
replace soz with ssz 
endif 
fx=sox 
fy=soy 
f z=soz 
skip 
loop 
else 
exit 
endif 
enddo 
exit 
else 

skip -1 
loop 
endif 
endif 
enddo 

skip -1 

ix=sox 

do while ,t. 

if sox= (ix-0 . 05) .and. (recno()ol) 
skip 
iy=soy 
skip 

do while .t. 

if iy<0.0001 

if soy>(iy+0.5) 
fx=sox 
fy=soy 
f z=soz 
skip -1 
do while .t. 

if soxo ( ix-0 . 05) 
sx=sox 
sy=soy 
sz=soz 

m= (sz-fz)/ (sy-fy) 
if ia>l 

ssz=f z+ (sy-fy) 
replace soz with ssz 
endif 
fx=sox 
fy=soy 
f z=soz 
skip -1 
loop 


else 

exit 

endif 

enddo 

ix=sox 

exit 

else 

skip 

loop 

endif 

endif 

if iy>-0 . 0001 

if soy<(iy-0.5) 
fx=sox 
fy=soy 
f z=soz 
skip -1 
do while .t. 

if soxo (ix-0 . 05) 
sx=sox 
sy=soy 
sz=soz 

m=(sz-fz)/ (sy-fy) 
if m<(-l) 

ssz=fz- (sy-fy) 
replace soz with ssz 
endif 
fx=sox 
fy=soy 
f z=soz 
skip -1 
loop 
else 
exit 
endif 
enddo 
ix=sox 
exit 
else 
skip 
loop 
endif 
endif 
enddo 
loop 
else 

if (recno()<>l) 
skip -1 
loop 
else 

iy=soy 

do while . t. 

if soy<(iy-0.5) 
fx=sox 
fy=soy 
f z=soz 
skip -1 
do while . t. 

if (recno()<>l) 
sx=sox 



sy=soy 

sz=soz 

m=(sz-fz)/ (sy-fy) 
if m<(-l) * 

ssz=f z- (sy-fy) 
replace soz with ssz 
loop 
endif 
fx=sox 
fy=soy 
f z=soz 
skip -l 
loop 
else 
exit 
endif 
enddo 
sx=sox 


sy=soy 

sz=soz 

m=(sz-fz)/ (sy-fy) 
if m<(-l) 

ssz=f z- (sy-fy) 
replace soz with ssz 
loop 
endif 
return 
else 


skip 

loop 

endif 

enddo 

endif 

endif 

enddo 

endif 

endif 

enddo 

return 



******************* PROGRAM NAME : DINFER1 . PRG ************************** 
USAGE : TO RAISE THE CUTTER TO A PROPER HEIGHT IF NECESSARY TO SOLVE 
* _ THE MACHING INTERFERENCE PROBLEM. 

^• PROGRAMMER : PAUL N. LIN DATE : DEC. 6, ‘1989 

****************************************************************** 
close databases 

select 1 


use dsoleto7 
select 2 


use dsoleto6 
ix=sox 


iy=soy 

iz=soz 

skip 

do while .t. 

if sox=ix.and. ( .not.eof () ) 
ux=sox 


uy=soy 

uz=soz 

rr=(3/32) 

pl=iz+rr 

p2=uy-iy 

p3=(p2*p2)+(pi*pi) -(rr*rr) 
if ( (4* (pl*pl) ) - ( 4 *p3 ) ) <-0.000001 
select 1 
append blank 

replace sox with ix,soy with iy,soz with iz 
select 2 
ix=sox 
iy=soy 
iz=soz 
skip 
loop 
endif 

zl=( (2*pl) +sqrt ( (4* (pl*pl) ) - (4*p3 ) ) ) /2 

if = zl>z2 1) ~ Sqrt( (4 *( pl *P 1 ) )"(4*p3) ) ) /2 

ZZ=Z2 

else 


zz=zl 
endif 
if zz<uz 

sigma=abs ( zz-uz) 
select 1 
append blank 

selict S 2 S ° X With ix,SOy with iy' soz wi th (iz+sigma) 
ix=sox 
iy=soy 
iz=soz 
skip 
loop 
else 

select 1 
append blank 

replace sox with ix,soy with iy,soz with iz 

select 2 

ix=sox 

iy=soy 

iz=soz 


skip 

loop 

endif 

else 

if soxoix.and. (.not.eof () 1 * 

select 1 w ' 

append blank 

select e 2 S ° X With ix,soy with iy.soz with iz 

ix=sox 

iy=soy 

iz=soz 

skip 

loop 

else 

select 1 
append blank 

replace sox with ix,soy with iy,soz with iz 

endif 

endif 

enddo 

close databases 
return 


PROGRAM NAME : DNCT 07 . PRG 




******************* 
close databases 
RE= ' $ ' 

STA= ’ (DNCT07 ) ' +RE 
select 2 
use dncto7 
APPEND BLANK 
REPLACE STAMENT WITH STA 
select 1 
use dsoleto7 
N='N' 

LI=100 
G1='G00 ' 

G2= ' G01 ' 

F='F60 ' 

; XX='X' 

YY= ' Y • 

ZZ='Z' 

PE=» (E) ' 

M1='M00 ' 

M2='M03 ' 

S 1=LTRIM ( STR ( LI , 4 ) ) 

S2=LTRIM ( STR ( DS0LET07->S0X ,8,4)) 
S3=LTRIM(STR(DS0LET07->S0Y, 8,4) ) 

S4=LTRIM (STR (DS0LET07->S0Z ,8,4) ) 

UP= ' 5 ' ' ‘ ‘ 

IF LEN (SI) <4 
S1=«0'+S1 
ENDIF 

STA=N+S1+PE+ ' G90 ' +RE 
SELECT 2 

APPEND BLANK 4 

REPLACE STAMENT WITH STA 

LI=LI+1 

se^ctT +LTEIM ' STR(LI ’ 4)) + ' (9) ,+ ' K06, + " r01,+RE 

, APPEND BLANK 
REPLACE STAMENT WITH STA 
LI=LI+1 

SELECT ' 2 ' +L ^ RIM (L ^ ,4) ) + l (9 ) , + l M°3 1 + ’ S7 57 1 +RE 

APPEND BLANK 

REPLACE STAMENT WITH STA 
LI=LI+1 

SELECT ' 2 ' +LTRIM (STR < LI - 4 > > + EE+G1+XX+S2 + VY+S3+Z2+S4 + RE 
APPEND BLANK 

REPLACE STAMENT WITH STA 

SELECT 1 

SKIP 

LI=LI+1 

S 1=LTRIM ( STR ( LI , 4 ) ) 

S2=LTRIM (STR ( DS0LET07 ->SOX , 8 4)) 

S3=LTRIM(STR(DS0LET07->S0Y, 8 %) ) 

S4=LTRIM(STR(DS0LET07->S0Z ,84)) 

IF LEN (SI) <4 'll 

S1='Q'+S1 

ENDIF 

STA=N+S1+PE+G2+XX+ S2+YY+S3+ZZ+S4+F+RE 

SELECT 2 


append blank 

replace stament with sta 

SELECT 1 
SKIP 

DO WHILE .T. 

IF . NOT. EOF () 

LI=LI+1 

S1=LTRIM (STR (LI , 4 ) ) 

S2=LTRIM (STR ( DS0LET07->S0X / 8,4) ) 
S3=LTRIM (STR ( DS0LET07->S0Y , 8 4)) 
S4=LTRIM(STR(DS0LET07->S0Z 8 4) 

IF LEN (SI) <4 ' 

S1='0'+S1 
ENDIF 

STA— N+S1+PE+G2+XX+S2+YY+S3+ZZ+S4+F+RE 

SELECT 2 
APPEND BLANK 

REPLACE STAMENT WITH STA 
SELECT 1 
SKIP 
LOOP 
ELSE 


LI=LI+1 

STA=N+LTRIM(STR(LI, 4) ) +PE+ » G01 ' +ZZ+UP+RE 
SELECT 2 
APPEND BLANK 

REPLACE STAMENT WITH STA 
LI=LI+1 


STA=N+LTRIM ( STR ( LI , 4 ) ) +PE+ 'GOO 
APPEND BLANK 

REPLACE STAMENT WITH STA 
LI=LI+1 


+XX+ '0.0 ' +YY+ '0.0 ' +ZZ+ '0.0' +RE 


STA-N+LTRIM ( STR ( LI , 4 ) ) +PE+ ' G17 ' +RE 
APPEND BLANK 

REPLACE STAMENT WITH STA 

li=li+i 


fnnJ^ LTRIM ^ STR ( LI ' 4 ) ) + ' (9) 1 + ' M05 ' +RE 
APPEND BLANK 


REPLACE STAMENT WITH STA 
LI=LI+1 


™ L Bsr R(Li ' 4),+,(9) ' + ' Mo2,+RE 

REPLACE STAMENT WITH STA 
LI=LI+1 


STA=N+LTRIM ( STR ( LI , 4 ) ) + 1 
APPEND BLANK 

REPLACE STAMENT WITH STA 
STA= 1 END » +RE 
APPEND BLANK 


(9) 


REPLACE STAMENT WITH STA 
EXIT 
ENDIF 
ENDDO 


1 + 'M3 0 ' +RE 


CLOSE DATABASES 


APPENDIX 3 


/* 

* Program s detachl .c 

* Purpose : Mainly to replace “detach. c" 

* Detach points below maximum width for each lattitude and reqist 

* them. 

* Date s 1 1 —Feb— 1990 

* Time : 14:35:42 

* 

* To detach: 

* Tor each lattitude find a point with smallest 
set maxX and minX 

* for each lattitude 

* step up from that point to both direction 

and register the vertices, 
update vertices and extreme values, 
for each lattitude 

write only vertices between the extrema. 

cal 1 ing : 

wsole<); writes sole array and initialize it. 
pbottom ( , ) ; 

wso lexyz ( , , ) ; 

*/ 

#include <stdio.h> 

^include <malloc.h> 
ttinclude <math.h> 

#def ine NUMLAT 59 
#def ine NUMLON 64 
struct coor C 
long x; 
long y; 
long z; 

struct lat C 

struct coor xy z CNUMLAT 1 ; 

1 5 

struct lat sit NUMLON 1 ; 
struct coor so 1 e C NUMLON 3 ; 

FILE *filein, tfileoutl, *fileout2; /* *fileout3, »fileaut<t; */ 
int lg, it; 

int xplus? xminusj lonp lus T lonminus; 


z 


to find maxX and minX 


main ( ) 

C 

printf (“enter read i n < ) \n“ ) ; 
read in( ) ; 

printf (“but read i n < ) and enter solext\n") ; 
or igen ( ) ; 
so lex t ( ) ; 

printf ("out so 1 ex t ( ) \n“ ) ; 

/^fcioseal 1 (); */ 



< 

filein = fopen C'ronXxX.rec", « r M ); 
for (lg = 0; lg < NUMLON; lg + +) < 

for (It = O; It < NUMLAT ; lt++) < 

fscanf (filein, "*/.D %D */D M , 

*-slClg3.xyzClt3.x, «.sl C lg 3 . xyzC 1 1 3 . y , 8.sl C lg 3 . xy z C 1 1 3 . z > 
3 

fc lose ( f i lein ) ; 

> 


or igen ( ) 

< 


/♦ 


/* 


/* 


/* 

> 


fileouta = fopen < " long . 1 at M , " w " ) ; 
fileout3 = fopen (“long.lip", M w”); 
for (It = 0; It < NUMLAT; lt++) ( 

for (lg = O; lg < NUMLON; lg++) C 

fprintf (fileouta, "Xld Xld Xld\n n , 
si Clgl . xyz C 1 1 3 .x, si C lg 3 . xyzC 1 1 3 .y , 
fprintf (fileouta, "Xld\nXld\nXld\n 
si Clg 3. xyzC It 3.x, s 1 C lg 1 . xyzC 1 1 3 -y , 


for lisp program */ 


slClg3.xyzClt3.z) ; 

3 

sl C lg3 - xyzC 1 t3 .z > ; 


fprintf ( f i 1 eou t3 , ,, n\n , *>; */ 

> 


f c lose (fileouta); 
fclose (f i leout3 ) ; */ 


*/ 


#def ine MIN ( long ) ( -3147^836^7 ) 

long zmax; 
struct cyxyz C 
long x; 
long y; 
long z; 

struct cyxyz *ptrnext; 
struct cyxyz *ptrprvs; 

> ; 

struct cyxyz *ptrfirst, *ptrthis, *ptrnew, »ptrZmax, *lring, *rrinq; 

so lex t < ) 

C 

printf ("start of so 1 ex t ( ) \n" ) ; 
zmax = MIN; 

fileoutl = fopen ( "so le. da t " , "w" ) ; 

/* fileouta = fopen < “sole. 1 ip” , “w” > ; // for lisp program */ 

xplus = 0; xminus = 0; 
lonplus = lonminus; 

printf ("solext loop entering\n" ) ; 
for <lt = 0; It < NUMLAT; lt++) C 
printf ("enter maker i ng () \n ") ; 
maker i ng ( ) ; 

printf ("enter w i deso 1 e ( ) \n" ) ; 
wideso le( ) ; 

printf ("enter edgead just < ) \n" ) ; 
edgead just ( ) ; 

printf ("enter wr i teso 1 e < ) \n ” ) ; */ 

wr i teso le ( ) ; 

/* ^ relmemO; //PENDING PENDING PENDING PENDING PENDING PENDING PENDING « 

) 


makering < ) 
C 



tpnntt < 1 x ieout 1 , M ma I lac error at It 7.d Xg %d\n a '»ig> it); 
ex i t ( 1 ) ; 

> 

zmax — MIN; 

ptrfirst = ptrthis = ptrnew; 
ptrthis->x = si CO 3 . xyz 1 1 1 3 . x ; 
ptrthis->y = s 1 CO 3 . xy z C 1 1 3 - y ; 
p t r th i s — > z = si C03 . xyzC 1 1 3 . z ; 
p trthis->p trnex t = {struct cyxyz *)NULL; 
ptr th is->p trprvs = (struct cyxyz *)NULL; 
ptrZmax = ptrthis; 

/* if (zmax < ptrthis->z) C 
// zmax = ptrthis->z; 

// ptrZmax = ptrthis; 

// 3*/ 

for (lg = l; lg < NUMLON; lg++) C 

if ( (ptrnew = (struct cyxyz *-) malloc (sizeof (struct cyxyz))) == NL 
C 

fprintf ( f i Ieout 1 , M mal loc error at It Xd lg 7.d\n",lg, It); 
ex i t ( 1 ) ; 

> 

p tr th i s->p trnex t = ptrnew; 

p tr new->ptrprvs = ptrthis; 

ptrthis = ptrnew; 

ptr th i s-> x = s 1 C 1 g 3 . xyz C 1 1 3 . x ; 

p tr th i s->y = s 1 C 1 g 3 . xy z C 1 1 3 . y ; 

p tr th i s->z = slClg3.xyzClt3.z; 

p tr th i s->p trnex t = (struct cyxyz *>NULL; 

/* if (zmax < ptrthis->z) C 

// zmax = ptrthis->z; 

// ptrZmax = ptrthis; 

// }*/ 

} 

ptrthis->p trnext = ptrfirst; 
p trf i rst->p trprvs = ptrthis; 

> 

wideso le< ) 

< 

long xthisi xthat, xmax, diff; 
struct cyxyz *this, *that; 

this = that = lring = rring = (struct cyxyz *)NULL; 

this = p tr Zmax— >p trnex t ; 
that = p tr Zmax— >p trprvs ; 
xthis = this->x; 
xthat = that— >x ; 
xmax = 0; 

/ * diff = labs (xthat - xthis); */ 
diff = abs (xthat — xthis); 

while (xmax < diff) < 
xmax = diff; 
lring = that; 
rring = this; 
that = that->p trprvs ; 
this = this->p trnext ; 
xthat = that— >x; 
xthis = this->x; 

/* diff = labs (xthat - xthis); */ 
diff = abs (xthat - xthis); 

) 

3 



«aetine HMAUtbHtt du.u 

^define RTOD 57,29578 

^define MAXRAD 0.3-4906585 

#def i ne ANGLE ( x , y ) < atan2 ( y , x ) ) 

#define SL0PE(x,y) ( (float) (y/x)) 

ttdefine PI 3.1415927 

^define M 0.363970 


edgead just ( ) 

< 

struct cyxyz *p0, *p 1 , *pm; 

long xO, xl» zO? z 1 , xm, zm, zadj, xwatchl, xwatchr ; 
float diffx, diffz, angle, slope; 

printf ("start of edgead just ( )\n M ) ; 

pO = pi = (struct cyxyz *)NULL; 


/* 


/* 

// 

// 

// 

// 


printf ("left edge\n") ; 

/* left edge */ 
pO = Iring; 
pi = Ir ing->p trnex t ; 
do C 

xO = pO— > x ; x 1 = p l->x ; 
zO - pO~> z ; z 1 = pi->z ; 
pO = pi; pi = pl->ptrnext; 

diffx = fabs ((float) (xl-xO)); diffz = fabs ((float) (zl 
if (diffx 1= 0) C 

angle = (float) atan2 ( ( doub 1 e )d i f f z , ( doub 1 e ) d i f f x 

> 

else 

angle = MAXRAD+1 ; 

> while (angle > MAXRAD); 
iring = pO; */ 

xwatchl = xO; 
da C 

pi = pO; 

pO = pO->ptrprvs; pm - pO->ptrprvs; 
xl = p 1 -> x ; xO = p0~> x ; xm = pm— > x ; 
zl = pl->z; zO = po-> z ; 
if ( xm <= xwatchl) ( 
xwa tch 1 = xm ; 

diffx = (float) (xl-xO); diffz = (float) (zl-zO); 
if (diffx [= 0) C 

slope = diffx / diffz; 

zadj = ( long) (slope * ( xm — xO) + zO); 

> el se 

zadj = zO; */ 

zadj = ( long ) (M * (xm— xO)) + zO; 
pm->z = zad j ; 

} else ( 

Iring = pO ; 
pm = 1 r i ng ; 

> 

> while (pm 1= Iring); 


printf ("right edge\n" ) ; 

/* right edge *-/ 
pO = rring; 
pi = rr i ng->p trprvs ; 
do C 

xO = pO->x ; xl - p 1 - > x ; 
zO = pO->z; zl = pl->z; 
pO = pi; pi = p 1 — >p trprvs ; 

diffx = fabs ((float) (xl-xO)); diffz = fabs ((float) 


-zO) ) 
> ; 


(zl-zO) ) 



/* 


/* 

// 

// 

// 

-// 


angle = ttlaat) atand ( (doudie)dittz, (doubie)dittx > 

else 

angle = MAXRAD+1 ; 

> while (angle > MAXRAD ) ; 
rring = pO; */ 
xwatchr = xO; 
do { 

pi = pO; 

pO = pO->ptrnext; pm = pO->ptrnext; 
xl = p 1 -> x ; xO = p0-> x ; xm = pm— > x ; 

2 1 = pl->z; zO = p0-> 2 ; 
if (xm >= xwatchr) < 
xwatchr = xm; 

diffx = (float) (xl-xO); diff 2 = (float) (zl-zO); 
if (diffx 0) C 

slope = diffx / diffz; 

zadj = ( long ) (slope * < xm — xO ) + zO) ; 

> else 

zadj = zO; */ 

zadj = ( long ) (-M * (xm-xO) ) + zO; 

pm— >z = zad j ; 

> else C 

rring = pO; 
pm = rring; 

) 

} while (pm ! = rring); 


} 


arc tan ( x , z ) 
long x, z; 

< 

if (x -= 0.0) C 

if (z >= 0.0) C 
return (PI); 

> 

else C 

return ( -P I ) ; 

> 

> else return (atan2( z,x )); 

} 


writesole( ) /*- with zigzag in mind! */ 

static int zigzag = O; 

switch (zigzag) <. 
case 0 : 

ptrthis = lring; 

while (ptrthis (= rring) i 

fprintf (fileoutl, ,,# 41d V * Id ttldNn", 
ptrthis->x, ptrthis->y, ptrthis->z); 
ptrthis = p tr th is->ptrnext ; 

fprintf (fileoutl, "%ld '/.Id */.ld\ n", 
ptrthis->x, ptrthis— >y, ptrthis->z); 



case 1: 

ptrthis = rring; 

whi le . ( ptrthis 1= Iring) < 

fprintf ( fi leout 1 * "7.1d 7.1d 7.1d\n'*, 

p tr th is->x > ptrthis->y, ptrthis->z) 
ptrthis = ptrthis->ptrprvs; 

> 

fprintf < f i leout 1 , "V.ld %ld 7.1d\ ** 
ptrthis— >x, ptrthis->y> ptrthis- 
break ; 

> 

if (zigzag == 0) < 

zigzag = 1 ; 

> else C 

zigzag — 0; 

> 


relmem( ) 

C 

int i ; 

struct cyxyz *this, *that; 


this = that = (struct cyxyz *>NULL; 

this = ptrfirst; 
that = th i s->p trprvs ; 
this = that->p trprvs; 
for (i = 0; i < NUMLON; i++) < 
free ((void *> that); 
that = th i s->p trprvs ; 
this = that— >p trprvs ; 

> 

> 


C A 


/* 

Program : insole. c 

Purpose : to find the inner sole by elevating SOLE . DAT and get difference 
between the data file and a plane. 

Date : 20-Feb-1990 
Time : 11: ^7 : 55 

*/ 

^include <stdio.h> 

^include <math.h> 

#include <malloc.h> 

struct cyxyz ( 
long x; 
long y; 
long z; 

struct cyxyz *next; 
struct cyxyz *prvs; 

>; 


struct cyxyz *head * *tail* *maxl, *max2; ; 
int numl inks = 0, 
numlat = 0; 
float z 1 , z2; 


ma in ( ) 

< 


*/ 


read in( ) ; /* 

z 1 z2 ( ) ; /* 

p 1 anequ ( ) ; /* 

path ( ) ; /* 


read in SOLE • DAT and save it into a linked list. */ 
find 2 maximum Z by dividing area of searching into 2 */ 
generate a plane equation with parameters provided by z 1 zl 
write a REAL sole derived by the above processes into a f 


> 


/* Function name : readin<) 

/* Purpose : Read SOLE.DAT in and make a iinked list whose structure is defi, 
3S the above global part. 

/**********•***********«*« a.***********/ 

read i n ( ) 

< 

struct cyxyz #now; 

long tmpx, tmpy, tmpz, refy; 

P 1 LE *f i 1 e i n ; 


filein = fopen ("sole.dat' 


’) ; 


tail _ , n ! W = <struct C Y*Y Z *> malloc (sizeof (struct cyxyz)); 
head->prvs = (struct cyxyz *) NULL; Y Y 

tail->next = (struct cyxyz *) NULL; 

fscanf (filein, '"/.D Y.D Y.D", &(head->x), 8.(head->y>, i(head-)z)); 
rety = head->y; 

numl inks++; numiat++; 

while (fscanf (filein, "%D */. D Y.D", *.tmpx, itmpy, itmpz) .= EOF) { 
now = (struct cyxyz *> malloc (sizeof (struct cyxyz)); 
if (now =— NULL) t 

printf ( error in now allocation. Exiting to system \n" ) ’ 
exi t ( 1 ) ; * 


now->prvs = tail; 
tail = now; 

tai l->next = (struct cyxyz *) NULL; 
tai l->x = tmpx; tail— >y = tmpy; tail->z = tmpz; 
if (refy 1= tmpy) C 
ref y = tmpy; 
num 1 a t++ ; 

num 1 i nks++ 5 

fclose (filein); 

> 

/* Funtion name : zlz2( ) 

/* Purpose : Search for 2 maximum z in two areas; first half and second half 
the bottom part of the foot. 

These 2 values are global and used in planequC) function. 

#define MIN ( long )( -2147483647 ) 
zlz2< ) 

< 

int count = 0, 
half lat ; 

long zmaxl = MIN, 
zmax2 = MIN, 
refy; 

struct cyxyz *naw; 

maxi = max2 = (struct cyxyz *) NULL; 
now = head ; 
refy = now— > y ; 
halflat = numlat / 2; 
while (count < halflat) <! 
if (zmaxl < now->z) ( 
zmaxl = now->z; 
maxi — now; 

> 

if (refy ! = now— >y) { 

refy = now-) y ; 
count++ ; 

> 

now = now— >next; 

/* printf ( "lowcount = 7.d\n M , count); */ 

> 

/* printf ("zmaxl %ld zmax2 7.1d\n" , zmaxl, zmax2) ; */ 

while ((count O numlat) && (now NULL)) i 
if (zmax2 < now->z) C 

/* printf ("zmax2 = 71d\n", zmax2); */ 

zmax2 — now-) z ; 
max2 = now; 

) 

if (refy !— now->y) { 
refy = now-) y ; 
count++ ; 

> 

now = now— >next; 

/* printf ("hi count = 7.d zmaxl 7.1d zmax2 %ld\n", count, zmaxl, zmax2); 

*/ 

/* printf ("zmaxl 7.1d zmax2 7.1d\n" , zmaxl, zmax2) ; */ 

> 

/* Function name : planequ() */ 

/* Purpose ; Generate a plane equation whose slope is 

slope = ( z2 + 25400) - (zl + 12700) 



Calculate zpiane value and then find difference with a point o 
the hee 1 - 

Replace the z of that paint by this difference. 

/**********♦**************-*-**--* ******** *■*■**♦ -*.*♦**■* «*«**«**•»•»* / 

p 1 anequ ( ) 

< 

int count; 

double zpiane? d i f f ; 
double slope; 
struct cyxyz *now; 
long y 1 , zl; 

slope = ( (doub le) ( (max2->z + 13700) - (maxl->z + 25400))) / 

( (double) ( (max2->y ) - ( max 1 ->y ) ) ) ; 

yl = max 1 — > y ; 
z 1 = max 1 ->z ; 

pr intf ( “s 1 ope = %f \n" , slope); 

now = head ; 

while (now •* NULL) < 

zpiane = slope * ( (double) (now->y ) - (doub le ) ( y 1 ) ) 

+ ( (double) < zl ) + (double) (25400) ) ; 
printf (“slope = %2-4f zl,yl %ld XI d zpl = Xld now->y = Xld now->z 

i 

slope ,max 1— >z ,maxl— >y, ( long ) zp lane, now— >y, now->z) ; 
diff = zpiane — (double) (now->z) ; 
printf (“diffz = XldXn", (long)diff); 
now->z = (long) diff; 
now = now->next; 

> 

> 


path ( ) 

C 

struct cyxyz *now; 

FILE *f i leout ; 

fileout = fopen ( "path.dat" , “w" ) ; 

/+ f printf (fileout, *'Yld y.d\n“ , num links, numlat) ; */ 
now = head; 
do t 

/* if ((now == maxi) !! (now == max2) ) C 

// fprintf (fileout, "Xld */.ld */• Id max\n“, now— >x, now->y, now— >z); 

// > else C */ 

fprintf (fileout, “Xld Xld '/.ld\n“ , now— >x, now->y, now->z); 

/* > */ 

) wiuie i now •— null;; 


> 




